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THE C HEMICAL BINDING EFFECT ON lOIUTRON THEHMALIZATION. 

INTRODUCTIOH 

Computlijg slow neutron spectra in reactors one must take into 
account that slow neutron differential cross section are essential- 
ly influenced, by the interatomic chemical bonds. This problem was 
solved in the past by means of simplified models. How, owing to the 
development of slow neutron spectroscopy, it became possible to ob- 
tain from experiments necessary information on bound atom dynamics 
and, making use of previously developed theoretical methods,^ to com- 
pute adequately slow neutron differential cross- sec tions.lt is well 
known, that a fine structure of differential cross-sections has only 
a minor influence on slow neutron spectra, so to calculate these it 
is quite sufficient to use an approximation which does not pretend 
to describe the fine structure but gives correct group constants with 
any reasonable choice of group widths. Such approximation, universal- 
ly adopted nowadays, is the incoherent gaussian approximation, which 
assumes that neutron waves scattered by different atoms do not inteiv 
fere and the space-time self-correlation function of an atom has a 
gaussian form. 

In this report the methods of, and computer programs for, calcu- 
lation of slow neutron group constants now used for reactor studies 
are described. The results of calculations for some moderators are 
presented and the influence of chemical binding on the integral pa- 
rameters important for neutron thermalizatiom, is studied. To analy- 
ze the influence of chemical bonds on neutron thermal izat ion neutron 
spectra in infinite media are studied, as well as transient phenome- 
na in space and time. The theoretical results are compared with the 
experimental data what may serve as a basis for the justification of 
one or an other model of the atom dynam ics, 

^Systematic treatment of the theory of slow neutron scattering 
together with the bibliography may be found in / 1/, 
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PROGRAM "UPRAS" 

Por multigroup differential cross section calculations a corn- 
put er program UPRAS ("Universal Program for Cross Sections Calcula- 
tions") was written in AKtOL- 60 language and autho mat ic ally trans- 
lated into one of the machine languages. 

On "bhe first stage^f its work this program calculates the sym- 
metric scattering law S(c^,t3 , three modes of atomic motion being 
taken into account: 

1, Acoustic vibrations with arbitrary spectinim. For small 9^ the 
phonon expansion is used^ f or large the Egel staff and Schofield 
formulae /2/ obtained by the steepest descent method may be used, 

2, Optical vibrations with a zero width spectrum, (Optical vibra- 
tions with a spectrum of a finite width can be treated as a part of 
acoustic vibrations), 

3, Diffusive motion with a half-width of the gaussian self-corre- 
lation function growing proportionally |t|' for large |t]'. In the ge- 
neral case the scattering law is computed as a convolution of the 
scattering laws for specified modes with arbitraiy weight coeffici- 
ents, The scattering law calculated, the matrices 
of zero and first angular moments of the double differential cross- 
sections are computed by integration. The matrices may be up to 
50x50 order. Integral chara cteristics of the differential cross-sec- 
tion, such as 'Oj , ^ , cos 0 eto,,also are calculated. 

The total cross section and thermal neutron sources may 

be as well calculated by the program UPRAS to give full infomiation 
for slow neutron spectra oaloulatlons, 

PROQHAM "PRASSIV" 

For the comparison of theoretical predictions with experimen- 
tals results on the scattering of slow neutrons in solids and liqui- 
ds a program PRASSIV ("Program for Calculation of Cross Sections by 
Time Integration") was written, wich calculates in the incoherent 
gaussian approximation the double differential cross-sections of 
slow neutron scattering by a method, entirely different from the one 
described above. 

In this program the symmetric scattering law is oalonlated by 
direct time integration of the well known Fourier-transformation 

36 0 

- 2 - 


Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 0011 0008-1 





Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 0011 0008-1 


expression with a gaussian cut-off factor introduced 

into tl^e integrand, '^e scattering law so calculated 

IS, due to a property of j;ourier transformations, the convolution of 
the true pattering law with the gaussian .of the width 

but this IS just what one measures experimentally! The quantity B 
determined by the resolution of the experimental device is speci- 
fied in the program as a function of and S ( for B as a matter 
of fact strongly depends on and E).So the presence of the cut- 
off factor not only enables to calculate the scattering law by sim- 
ple time integration, but also takes into account the experimental 
resolution when comparing the theory with experiment. 

If there is no diffusive motion the function ^ Ct) is calcu- 
lated through the normal frequency spectrum . The spectrum is 

arbitrary and the presence of optical vibrations does not compli- 
cafe calculations. If the diffusive motion is present the function 
is calculated by the interpolation formula proposed previous- 
ly /V» (at small t the diffusion motion does not essentially change 
the function in this model). 

The program PRASSIV can compute double differential cross-sec- 
tions G'(,Eo~~*E, 03 as well as differential cross-sections (TtEp-* E) 
and and their integral characteristics. It can be used 

for multigroup matrices calculations if only B is taken not larger 
than the group width. 

Differential cross-sections. 

We shall briefly characterize the problems arising in the cal- 
culation of differential cross-sections for some important modera- 
tors and present some results of calculations, 

frequency spectrum of the vibrations of a hydrogen atom 
in a water molecule can be represented as a sum of two terms. 

The first ("acoustic") one di scribes acoustic vibrations and hinde- 
red rotations of the molecule and extends up to approximately 0,1 3ev, 
The second ("optical") term describes the vibrations of the hydro- 
gen atom in the molecule, characteristic frequences being "^0,2 ev 
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and <^0,48 ev, Krieger and Nelkin /3/ found the weight of the 
"acoustic” term A =0«48» This value was used in our calculations. 

EL 

The form of the "acoustic" spectrum was determined by Larsson /4/ 
from the scattering of cold neutrons and by Egel staff /5/ and Mos- 
tovoy /6/ from double neutron spectrometer measurements by means of 
the Egelstaff extrapolation technics /7/. 

The differential cross-section S' with initial energy 
E^=O,102 ev was calculated by the program PRASSIV, two frequency 
spectra being used: those by Larsson and by Egelstaff. The results 
were compared with the experimental results of Mostovoy (see Pig.l). 
One can see that the cross-section calculated with the Egelstaff 
freguency spectrum agrees well with experiment, so this spectrum was 
adopted for differential cross-section calculations to be used for 
neutron thermalization studies in water. 

The integral characteristics of differential cross-sections are 
less sensitive to the foirm of the "acoustic" spectrum than differen- 
tial cross-sections themselves, wich can be seen from figures 2 and 3 
where the mean logarithmic energy loss \ and the average cosine of 
the scattering angle tos © are shown, Theoretical curves for both 
spectra agree well with experiment. 

Graphite Various theoretical approaches give different frequency 
spectra of graphite lattice vibrations. The experimental determina* 
tion of frequency spectrum from neutron scattering data by the Egel- 
staff extrapolation method is not very reliable because of strong 
coherence of the scattering. But for reactor studies a model which 
gives, accurate enough, the integral characteristics of scattering 
cross-sections is quite satisfactory. Such a model was proposed by 
Egelstaff /8/. Using this model, some investigations were carried 
out of neutron thermalization in space and time, which will be des- 
cribed below. Besides, the differential scattering cross-section was 
calculated with a frequency spectrum obtained by Baldock /9/. 

The integral characteristics of the differential cross-sections cal- 
culated by these two models were compared. The distinctions are not 
very large. The mean energy loss is shown in Pig. 4 as an example. 
Beryllium, The high symmetry of beryllium crysal lattice affords to 
use the Debye model with the Debye temperature 0 =1000°K found from 
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heat capacity measurements. In Pig. 5 the differential cross-section 
GCe.o-»E) for the Debye crystal with the atom mass 9 is shown in 
order to demonstrate the work of the program PPASSIV when energies 
are large. The temperature T was taken equal to the Debye tempera- 
ture 0 . The energy E is measured in initial energy xuiits. 

Por =» 60 0 the differential cross-section approaches to the pla- 
teau containing curve specific for the soattaring on a free atom. 
Zirconium hydride. The frequency spectrum of vibrations of a hydro- 
gen atom in zirconium hydride has an optical part with the average 
energy about 0,13 ev and the half-width about 0,025 ev and an acous- 
tic part with the Debye temperature about 0.02 ev /I0/,/537* 

The weights of these are aproximately in the ratio of the Zr mass to 
the H mass. This infomnation is sufficient to construct an approxi- 
mate frequency spectrum. ^ means of the program PEASSIY one can use 
experimental double differential cross-section (at one angle and ini- 
tial energy only) to correct a frequency spectium so as to fit the 
experimental points. The energy distribution of neutrons with initial 
energy 0,020 ev scattered at the angle 80° which was obtained by V.A, 
Parfenov et al. on the double neutron spectrometer using pulsed fast 
reactor in Dubna /11/, is shown in Pig. 6. 

The theoretical curve was computed by the program PEASSIV with the 
corrected frequency spectrum. 

lithium hydride. According to /23/ the frequency spectrum cf vibra - 
tions of a hydrogen atom in lithiiun hydride was taken as consisting 
of two optical lines; 0,10 ev and 0,13 ev with weight ratio 2 : 1 
and an acoustic part with the Debye temperature 9 =250°K and the 
weight My + Pig. 7 the total scattering cross-sec- 

tion in lithium hydride computed by the program "UPEAS” is shown in 
comparison with the experimental points of Doilnitsyn et al. /12/. 

Spectra in infinite media. 

The chemical bond effect on the neutron thermalization can be 
appreciated by the investigation of its influence on various inte- 
gral characteristics of neutron spectra. In particular, one of the 
most important thermalization characteristics of a moderator is the 
neutron gas temperature in an infinite media with not very large Vv 
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abeorption, when the enengy distribution ot thermal neutrons can 
be treated as a maxwelllan one. 

Por three hydrogenlous moderators: water, zlrconitun hydride and 
lithium hydride, the neubron spectra in an Infinite medium at 
To300°E poisoned by the Vv —absorber were computed for several val- 
ues of absorption cross-section varying up to 10 bam per hydrogen 
atoStThe spectra in water and llthiiun hydride are nearly marwellian 
with a certain effective "temperature T^*The spectra in zirconium hyd- 
ride (see Pig* 8) are essentially non-maxwell ian in the vicinity of 
vibrational quan'tum energy, but still it is possible to determine the 
effective temperature from the lower energy region. 

In order to represent in a convenient way the dependence of ef- 
fective neutron temperature on the absorption cross-section, we 
small make use of the formula 
Tk/t = 

where Oq, is the absoarption cross-section per hydrogen atom at the 
energy E a t, and a 20,4 b* is total scattering cross-section on 
the free hydrogen atom. The coefficient oC can be understood as an 
effective mass of moderator: oC. = Vy\ ,for in the moderation of neu- 
trons by an ideal monatomio gas is approximately equal to the mass 
of the atom /1 3/* The dependence of yy\^on for the three moderators 
is shown in Pig, 9, The strong dependence of YWe^on for 
zirconium hydride follows from the feet that the pattern of the slow- 
ing down process in zirconium hydride is quite different from that 
in a heavy gaseous moderator. On the contrary, the nearly constant 
auid close to unity value of Vv\ 4 ^ for water testifies that as it was 
noticed before by its' slowing down properties water is close to free 
hydrogen atoms. Por water and lithium hydride the calculated results 
agree well with the experimental data /12/,/14/. 

The rethermallzatlon length. 

In dealing with reactors with large temperature gradients it is 
important to know at what distance from the temperature discontinuity 
surface a neutron speotnun with Ideal j temperature is settled • 

This distance is characterized by so called thermalization or 
rethermallzatlon length which for the plane geometry is de- 

fined through the relation 
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T>^ (.OO'i = Ip VU tVv\ 

where the neutron gas temperature at the distanceVt^roa the 

temperature discontinuity plane. The rethermalization length is 
essentialy dependent on the chemical honds between moderator atoms 
and is one of the important characteristics to judge on the chemi- 
cal bond effect. As an example we present the results of a calcula- 
tion of neutron space— energy distribution in a graphite cylinder 
consisting of two zones with different temperatures. The inner zone 
with a radius of 40 cm has a temperature of 300 K, The outer zone 
with a radius 80cm has a temperature of 450®K. In Pig. 10 th® neutron 
gas temperature T^ against the distance from the cylinder axis 
is shown j calculated with the ideal gas model (M=42) differential 
cross-sections and with cross-sections computed by the program UPRAS 
with the Egel staff vibration frequency spectrum /sj. For free gra- 
phite atoms L\.'^=4,2 cm, fcr bound atoms Lt\^=6,9cm, As the retherma- 
lization length for ideal monatomic gas increases with M as fiA" we 
may say that the effective mass of a bound graphite atom is about 32. 
The rethermalization calculations were carried out by the program 
T-2 CF-1 approximation, 40 energy groups). 

Pulsed neutron therm al ization, 

A great deal of information on slow neutron scattering by bound 
atoms is contained in pulsed neutron source experiments. So the com- 
parison with these experiments is a good check for a theory. 

For this purpose a program '» TEPI" C Thermalization and Diffusion) 
was written in ALGOIr-60 language, wliloh calculatee in the P-1 approx- 
imation the time dependent energy distributions of neutrons in a 
finite moderator block. A system of two differential equations for 
neutron density IT(E,t) and flow j(E,t) is solved by the Runge— Kut— 
ta method for a space component with, given The maximum number of 
energy groups is 30. 

The time dependence of the average neutron velocity in water 
is shown in Fig.ll. In about 12 moseo after the pulse the effective 
neutron temperature T^^ begins to approach the media temperature T 
exponentially: T^^ -T = const » 

wi1;li a 'thenna.liza'tion ‘tiniB me S 80 # 

3 6 0 -7- 
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For graphite the exponential dependence settles at about 
800 no seoy the thermallsatlon time being equal to 315 me seo 
(at the density of 1,6 gr/on^), 

A special attention was paid to the investigation of asymp- 
totlo neutron spectra at large time when ir(E,t) n nC£)e f 
J(E,t) « l(E)e'‘<t. 

2 

The dependence of the decay constant ^ on B was measured 
for various moderators by many authors. From these measurments 
the diffusion ooeffiolent and diffusion cooling C were de- 
termined. It is well knoim that the equation for the asymptotic 
energy distribution of neutrons in the stationary diffusion pro- 
cess differs from that in the pulsed source case only by the si- 

2 

gn of buckling B • 

Measuring diffusion length as a function of absorption cross- 

-section it is possible to find out cC as a function of B in 
2 

the negative B region. Thus the asynptotlo speotra and the oon- 

stant((were calculated by the program TEBI both for positive and 

2 

negative values of B * 

The calculated asymptotic density and flowoeutron speotra sore 
shown in Pig, 12 for water and Pig, 13 for graphite at different 

p 

values of B . For both moderators the density and flow speotra 

differ appreciably. The effect of diffusion cooling for B >0 

2 

and diffusion heating for B < 0 is olearly seen. It may be wor- 
th-while to note at the energies larger than Tn(E) and l(E) are 

maxwelllan with some effective temperature T depending on B , 

2 ^ 

In Pig, 14 the function (B ) f6r water at the temperatxire of 
295®K is shown. The differential cross- sections were oaloulated 
by the program UPRAS with the frequency speotnua of Egelstaff 
The experimental points are taken from (l53andCl6j (pulsed 
source measurement^ and (j7j (poisoning experiments). The agreeme- 
nt of the theory with poisoning experiments (negative B ) is 

2 

very good. In the positive B region the theory agrees well with 
the measurements by Lopez and Beyster jj^and somewhat diverges 
from those by Antonov and Ktlchle (l^(Euohle*s date are not shown 
in the figure). 

We must note that die experiments cf lopsz and Better were osusied out 
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in a rectangular geometry while the experiments of Antonov and Ktt- 
ohle - in a cylindrical geometry in which case the determination of 
B is less soxind. The calculated values of the parameters are: 

1)^ =i 66 • 10^ , C = 3.0 S • 10 ^ 

On the whole one may conclude that the use of the Egelstaff frequ- 
ency spectrum for water leads to a good description of neutron 
thermalization and diffusion. 

The function for graphite at the temperature of 295°K Is shown 
^ ^^S»^5.The differential cross-sections were calculated ly the 
program UPRAS with the model frequency spectrum of Egelstaff /8/. 
Since there is no data on the poisoning method, the calculations 
were carried out only for 0. The experimental points are from 
/19/ and /20/.Por large dimensions of moderator bloohs the theory 
agrees well with the experiments, which Implies that the calculated 
diffusion coefficient 

is very close to the experimental value. Por smaller dimensions the 
theory gives larger , whloh implies that the theoretical value of 
the diffusion cooling constant C= is smaller than the 

experimental value. 

CQgCmsiOH 

Owing to the intensive work on the neutron thermalization prob- 
lem carried out in various countries the experimental and theoretloal 
methods are available at present, which make possible to obtain the 
differential cross-sections of slow neutron scattering on bound atoms 
and check them by comparison with experimental data on neutron speo- 
tra.Added to the up-to-date computational methods (see for instance 
/21/,see also /22/) this makes possible for any system of interest 
to calculate the space-energy distribution of slow neutrons ef- 
fects resulting from its variation. 
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Pig, I, The differential scattering cross->8ectlon In 

water at Bo* 0,I02ev»The solid curve is calculat- 
ed with the Egelstaff frequency spectruin/5/,the 
dotted one with the Larsson spectrum/4/. The ex- 
perimental points are from Mostovoy et al /6/. 



Pig. 2. The mean logarihlimic ene- 
rgy loss \ of neutrons 
in water. The solid curve 
is calculated with the 
Egelstaff spectrum, the 
dashed one with that of 
Larsson. The ei^erimental 
points are from /6/. 
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Fig. 3* The average cosine of 
the scattering angle of 
neutrons in water. The so- 
lid curve is calculated 
with the Egelstaff spectrum- 
The dashed curve with Lars- 
son's, the dot-dashed curve 
the Nelkin model. The expe- 
rimental points are from 6 
(empty circles) and 8 (fil- 
led circles). 



Fig. 4, The mean energy loss in graphite with the fre- 
cuency spectrum of Egelstaff 8 (solid curve) and 
Baldock 9 dashed curve). 
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Pig, 7* The total scattering cross-section on a hydro 
gen atom in lithimn hydride. 


• The netroun spectra in an infinite medium of 
zirconixmi hydride poisoned by i -absorber at the 
room temperature. The figures at the curves denote 
the absorption cross-section at E s 0.025 ov in 
bams per a hydrogen atom. 
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Fig. 9. The effective mass of -three hydrogenlous mode- 
rators, 0 - the experimental value for water 14 
6 - that for lithiium h^^ide 12 * 



Fig* IX>. The neutron gas tempena-ture in a graphite cy- 
linder consisting of -two zones with the tezo- 
pera-tures of JCX)® K AND 450® K» The solid curve 
is calculated with the bound atom differenti- 
al cross-section, the dashed curve by gas model 
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Pig. 12. The asymptotic neutron density and flow distribution 
in water 
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Pig, 14 1 Th.e decay constant as a ftmction of water 

X- experiment 1 17 1 ; o— experiment lI5]» 

• - experiment llSJ. 


pul (• 16 vJ 

T’lttr* 


Fig, 15. Tlie decay constant ac &s & function of B for graphite 
X- experiment 1191 5 O” experiment 120], 
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